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We calculate the temperature dependence of the transport properties of heavy-fermion systems such as resis- 
tivity, optical conductivity, thermoelectric power, the electronic part of the thermal conductivity, and the "figure 
of merit." The one-particle properties of the periodic Anderson model are obtained within dynamical mean-field 
theory for the paramagnetic phase using Wilson's numerical renormalization group and the modified pertur- 
bation theory as impurity solvers. We discuss the dependence of the transport properties on the band filling, 
valence, and Coulomb correlation U . The typical experimental findings can be reproduced and understood, in 
particular the temperature dependence of the resistance and the thermoelectric power and their absolute magni- 
tude for both metallic heavy-fermion systems and Kondo insulators. For large values of U , we find a negative 
Seebeck coefficient S{T) for an intermediate-temperature regime as observed in S{T) of CeCu2Si2. We ana- 
lyze different estimates for possible characteristic low-temperature scales of the lattice. Our results indicate a 
one-parameter scaling of thermodynamic and some transport properties with a strongly occupancy-dependent 
scaling function. This is consistent with a strong-coupling local Fermi-liquid fixed point of the effective site 
governing all low-lying excitations for T ^ in the paramagnetic phase. 

PACS numbers: 71.27.+a, 71.10.Fd, 72.10.-d, 72.15.-v 
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I. INTRODUCTION 

For nearly 30 years the investigation of heavy-fermion sys- 
tems (HFSs) has been one of the most fascinating and inter- 
esting fields in condensed matter physics both experimentally 
and theoretically.''^ The heavy-fermion phenomenon exists in 
a number of lanthanide and actinide compounds and mani- 
fests itself in the apparent existence of quasiparticles with very 
large effective mass m* at low temperatures T . This is seen al- 
ready in the specific heat, the electronic part of which shows 
the usual linear temperature dependence c = yr at very low 
T for many metallic heavy-fermion systems but with a y co- 
efficient being 100-1000 times larger than in usual metals. 
However, this y coefficient itself is strongly T dependent and 
rapidly decreases to "normal" values at higher T . Similarly, 
the static magnetic susceptibility x{T) crosses over from a 
Curie law xiT) ^ l/T at high temperatures to a Pauli-type be- 
havior for very low T: it approaches a finite value x{T — 0), 
which is also strongly enhanced compared to the Pauli sus- 
ceptibility of normal metals; the ratio between y and % (Wil- 
son ratio) is of order 1, indicating that both enhancements are 
caused by the same physical processes. 

The heavy-fermion behavior is driven by the local mo- 
ments of incompletely filled / shells of the lanthanide or ac- 
tinide ions. At high temperature, the weakly coupled / elec- 
trons cause mainly incoherent and with decreasing tempera- 
ture logarithmically growing spin-flip scattering for the con- 
duction electrons. Below a characteristic temperature scale, 
a crossover to a coherent low-temperature phase is observed. 
The / electrons contribute significantly to the formation of 
heavy quasiparticles while their moments are dynamically 
screened. Some heavy-fermion compounds such as SmBg 
or Ce3Bi4Pt3 exhibit similar behavior at high temperature, 
but for low T a crossover to an insulating heavy-fermion 



ground state with a narrow gap at the Fermi energy is found. ^ 
These heavy-fermion systems are termed Kondo insulators 
even though the insulating behavior is driven by Fermi vol- 
ume effects. In this case, the strong correlation induces only a 
narrowing of the band gap rather than the insulator transition 
itself as in Mott-Hubbard insulators. 

One particularly fascinating aspect of heavy-fermion com- 
pounds is that despite the existence of local moments they can 
undergo a superconducting phase transition, as observed'* in, 
e.g., CeCu2Si2. In other systems magnetic phase transitions 
are reported, with a large variety of types of order depending 
on the composition. ''^ In the weak-coupling or local-moment 
regime, the local / moments tend to order antiferromagneti- 
cally, while in the Kondo regime itinerant magnetism with an 
incommensurate ordering vector Q and strongly reduced or- 
dered moments is observed.^ This crossover from local to itin- 
erant magnetism has been studied in CeCu^_j)2Gs2 by vari- 
ation of the coupling constant upon doping and has been at- 
tributed to the strongly temperature- and coupling-dependent 
residual quasiparticle-quasiparticle interaction. ^'^'^'^ 

In this paper, we focus on the calculation of transport prop- 
erties of heavy-fermion systems in the paramagnetic phase. 
Characteristic information on heavy-fermion systems is ob- 
tained by measurements of the transport coefficients. The 
Seebeck coefficient S{T), or the thermoelectric power, is of- 
ten nonmonotonic, can exhibit different extrema, and even 
sign changes'"'"''^''^"'^''^ related to particle-hole asymme- 
tries. Much of the recent interest in the heavy-fermion ther- 
moelectricity is stirred by its large narrow peak at low temper- 
atures which might be useful for solid-state cooling devices. 

For high T, the resistivity p{T) is determined by a nega- 
tive temperature coefficient, and one usually observes a log- 
arithmic, "Kondo"-like increase of the resistivity p[T) with 
decreasing T. In Kondo insulators, p{T) crosses over from 
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low values to an activation behavior for smaller T, reflect- 
ing the insulating ground state. In metallic heavy-fermion 
systems — for example, CePdg (Ref. 17), CeAls (Ref. 18), or 
CeCug (Ref. 19) — a maximum is observed in p {T) before the 
resistivity approaches a small residual value for T ^ 0. At 
low temperatures T, p{T) often obeys a law in such mate- 
rials. The logarithmic increase of p{T) above a characteristic 
low-energy scale T\o^ (for a discussion of low-energy scales 
see Sec. VI) is related to growing spin-flip scattering of the 
conduction electrons with decreasing temperature as a mani- 
festation of the Kondo effect. Since the resistivity must vanish 
for r ^ in a translational-invariant system with a nondegen- 
erate ground state, a maximum in p{T) must connect these 
two regimes. The observation of a behavior well below 
the characteristic temperature scale Tjow in combination with 
a strongly enhanced y coefficient of the specific heat indicates 
the formation of a Landau Fermi-liquid by the heavy quasi- 
particles. 

We show in a comprehensive study that the low- 
temperature transport properties of heavy fermions can be un- 
derstood in terms of a minimalistic model of interaction be- 
tween local and itinerate degrees of freedom within the dy- 
namical mean-field theory. ^"'^^ Spin-flip scattering between 
these degrees of freedom yields a logarithmically increasing 
resistivity with decreasing temperature. Our calculations re- 
produce the experimentally observed maximum and the 
behavior at low temperatures characteristic of coherent trans- 
port in the Fermi-liquid phase. The calculated optical conduc- 
tivity exhibits a very narrow Drude peak and an optical ex- 
citation gap characteristic of heavy-fermion materials. The 
midinfrared peak is located in the correct frequency range but 
with a too narrow width compared to the experiment^^ due to 
the lack of local crystal electric field excitations in our model. 

We find very large absolute values for the thermoelectric 
power exceeding 150 /iV/K close to the Kondo insulator 
regime at temperatures of the order of Tiow • In this case, the 
purely electronic figure of merit can be larger than 1 . The ab- 
solute magnitude of the thermoelectric power agrees very well 
with the typical experimental reported values. 
The sign changes of ^(r) depend very sensitively on the 
particle-hole asymmetry and the band filling. Also in the ex- 
periment, the details of the thermoelectric power vary strongly 
between different materials. '"'"''^■'^■'^''^ 

The experimental evidence^'^^ compiled over the past ten 
years also indicates that even for heavy-fermion systems with 
a paramagnetic ground state, the temperature dependence of 
the specific heat and the magnetic susceptibilities often do 
not agree with the predictions of Fermi-liquid theory. There- 
fore, the phenomenological term "non-Fermi-liquid" was at- 
tributed to such regimes appearing in a large variety of differ- 
ent materials. ^'^^ Despite a tremendous experimental and the- 
oretical effort it is, however, still not clear whether the non- 
Fermi-liquid effects observed in heavy-fermion compounds 
are related to novel low-lying nonlocal excitations in con- 
centrated systems, true local non-Fermi-liquid physics^^'^^ or 
simply competing local energy scales. 

The physics in heavy-fermion compounds is driven by 
the interaction between two distinct subsystems: localized. 



strongly correlated / electrons hybridizing with extended con- 
duction bands. The periodic Anderson model (PAM) takes 
these ingredients into account, comprising of spin degener- 
ate conduction electrons, a lattice of correlated localized / 
electrons and a hybridization (cf. Sec. II A). The conduc- 
tion electrons experience spin-flip scattering from the mag- 
netic moments of the / shells. On the other hand, the model 
also accounts for the RKKY interaction between two localized 
magnetic moments in the particle-hole channel,^'^ leading to 
the competition between the Kondo screening and magnetic 
ordering.'^ While these qualitative aspects have been known 
for almost 30 years, a reliable calculation of the full phase di- 
agram of the PAM in more than one dimension has not been 
possible yet. 

In this work, we apply the dynamical mean-field theory^"'^' 
(DMFT) mapping the PAM on an effective single impurity 
Anderson model^^ (SIAM), which becomes exact^"-^' in the 
limit of infinite spatial dimension^^ for lattice models of cor- 
related electron systems. As impurity solvers for the effective 
SIAM we use Wilson's numerical renormalization group^^ 
(NRG) for all interaction strength and the modified perturba- 
tion treatment^"*'^^ (MPT) for small values of U. The NRG 
provides an accurate and non-perturbative description of the 
effective site^^ and yields the correct low-energy scale Tk for 
the SIAM (cf. Sec. VI B): Tk °^ exp(-l/g), with the dimen- 
sionless coupling g = 7po(0), the spin-flip scattering rate J, 
and po(0) being the density of states of the noninteracting 
conduction band of width D at its band center. Furthermore, 
the NRG provides detailed information on the possible fixed- 
point structure of the effective site. Compared to other possi- 
ble and frequently used DMFT impurity solvers such as quan- 
tum Monte Carlo^^ (QMC), exact diagonalization (ED), and 
noncrossing-approximation''^-^^'"^'^ (NCA), the NRG and MPT 
have the advantage to be applicable also for very low tem- 
peratures. This seems to be essential for the reproduction of 
Fermi-liquid behavior and simultaneously identifying the rel- 
evant low-temperature scale. 

The transport properties are calculated in the linear re- 
sponse regime using Kubo formulas. They relate the quanti- 
ties of interest — e.g., the frequency-dependent conductivity — 
to the charge and heat current-current response functions, 
which are two-particle Green functions. These response 
functions greatly simplify in any local approximation since 
current-operator vertex corrections vanish'*' even in the pres- 
ence of crystal field levels:^''^^ only the free particle-hole 
propagator enters. From the /- and band-electron self- 
energies calculated within the DMFT-NRG and DMFT-MPT, 
the static conductivity is obtained from the transport integrals 
in the limit o ^ 0, as well as the thermoelectric power S{T) 
and the electronic part of the thermal conductivity k{T). Re- 
cently, a very comprehensive study of the transport properties 
in heavy-fermion systems as well as Kondo insulators using 
the local-moment approach (LMA) has been published in a se- 
ries of papers by Logan and collaborators. Previously, 
Costi and Manini'*^ investigated the low-energy scales and 
temperature-dependent photoemission in the 5=1/2 Kondo 
lattice model using the DMFT-NRG. 

The scope of our paper is to calculate the resistivity, opti- 
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cal conductivity, the thermoelectric power, and the figure of 
merit for arbitrary interaction strength and band filHng within 
a single approach, the DMFT-NRG. We address the question 
of the existence of several low-temperature scales in the peri- 
odic Anderson model which might manifest themselves in the 
transport properties. The quality of the DMFT-NRG calcula- 
tions is critically examined by a comparison with the modified 
perturbation theory in the weak-coupling regime {U 0), ac- 
curately described by the DMFT-MPT. 

The paper is organized as follows: In Sec. II we introduce 
the model and the notations as well as the DMFT and the two 
impurity solvers used: the NRG and MPT. In order to gain 
physical insight into the temperature evolution of the trans- 
port properties, we discuss the single-particle properties such 
as the temperature-dependent /-electron spectral function and 
self-energy of the PAM in Sec. Ill, since they directly deter- 
mine the transport properties through the Kubo formulas. Sec- 
tion IV is devoted to the transport theory for the PAM and 
states explicitly equations for the different transport coeffi- 
cients investigated. We report results for the transport quanti- 
ties in Sec. V. The temperature dependences of the resistivity 
p(r), the thermoelectric power S{T), the thermal conductiv- 
ity k{T), and the temperature and frequency dependence of 
the optical conductivity a{co,T) are shown for different val- 
ues of U and band filling for metallic heavy-fermion systems 
as well as Kondo insulators. In Sec. VI we define precisely our 
different low-energy scales and discuss how they are linked. 
We conclude with a summary and outlook in Sec. VII. 



II. THEORY 
A. Model 

As mentioned above, the Hamiltonian of the simplest ver- 
sion of the periodic Anderson model is given by 

ka 

+ VY.iflcio+clfia). (1) 
ia 

Here, cj^ (cl ) destroys (creates) a conduction electron with 

spin C7, momentum^ and energy £j^. The energy Efa denotes 

the spin-dependent single particle /-level energy at lattice site 

f -j- 

;, fij^ = fiafi<^ /-electron occupation operator (per site 
and spin), /(j (/^) destroys (creates) an / electron with spin a 
at site /, and U denotes the on-site Coulomb repulsion between 
two / electrons on the same site ;. The uncorrected conduc- 
tion electrons hybridize locally with the / electrons via the 
matrix element V . While Eq. (1) includes possible Zeeman 
splitting of the energies in an external magnetic field H, we 
set // = throughout the remainder of the paper and treat all 
properties as spin degenerate. 

Even though only a single effective / level is consid- 
ered, this model is quite general. It describes any heavy- 
fermion system with odd ground-state filling of the / shell. 



for which in a strong crystal field environment the degener- 
ate Hund's-rule ground state may be reduced to an effective 
spin-degenerate Kramers's doublet. In addition, charge fluc- 
tuations to even / fillings leave the / shell in crystal field sin- 
glets. The Hamiltonian contains four energy scales. The inter- 
play between Ef and U controls the average / filling as well as 
the local moment formation for large U and negative £f. The 
Anderson width To = y^7rpo(0) determines the charge fluctu- 
ation scale of the / electrons with Po(0) being the density of 
states of the noninteracting conduction band of width D at its 
band center. 

The total filling per site, «tot = L<j((«;a) + ("fa))' is 
kept constant by a temperature-dependent chemical potential 
pL{T). We absorb the energy shifts into the band center of 
the conduction band, ej = + e^, as well as the / level ej. 
For viiot = 2 and U =0, the uncorrected system is an insu- 
lator at r = 0, since the lower of the two hybridized bands 
is completely filled. According to Luttinger's theorem a fi- 
nite U of arbitrary strength does not change the Fermi vol- 
ume which includes the full first Brillouin zone. As long as 
the ground state does not change symmetry due to a phase 
transition, the system remains an insulator at arbitrarily large 
Coulomb repulsion. Therefore, the nonmetallic ground state 
of Kondo insulators is not correlation induced, but it is already 
present for the noninteracting system and a consequence of 
Luttinger's theorem. For nonintegral values of «tot, the para- 
magnetic phase of the system must be metallic. 

In general, we can distinguish three different adiabatically 
connected regimes for U ^ Tq and £/ < 0. In the mixed 
valence regime, |e^|/ro ~ 1, the system is dominated by 
charge fluctuations yielding a nonintegral value of the / filling 
n / < 1 . In the stable moment regime | gy^ | ^ Fq , the / electrons 
remain strongly localized and form a stable local moment, 
which tends to order antiferromagnetically due to the RKKY 
interaction mediated by the conduction electrons. These two 
regimes are connected by the Kondo regime for moderate ra- 
tios je^l/Fo. The competition between screening of the local 
moment due to the Kondo effect and the RKKY interaction 
makes this crossover regime the most interesting one since it 
can lead to long-range magnetic order of the residual magnetic 
moments. Spin-density-wave, metallic and superconducting 
ground states are observed in heavy-fermion materials' which 
are believed to be described by the Kondo regime. 



B. Dynamical mean-field theory 

Setting aside exact solutions in one dimension^^ using the 
Bethe ansatz for the Kondo lattice model, to our knowledge 
no exact analytical solution has been found for the model (1) 
with finite U. Therefore, one has to rely on suitable approxi- 
mations for the PAM. An obvious first approximation is the as- 
sumption of a purely local, site-diagonal (i.e., ^-independent) 
self-energy, which for the PAM is even better justified than 
for other lattice models of correlated electron systems, as 
the first corrections are at least of order V^. Within a lo- 
cal self-energy approximation the complicated lattice prob- 
lem can be mapped on an effective impurity problem; i.e., the 



4 



PAM can be mapped on an effective SIAM.^^ Such a map- 
ping was first used already about 20 years ago in connection 
with appHcations of the NCA (Refs. 38,39,40 and 48) to the 
PAM7''*^'^"'^' -^^ a self-consistency condition accounts for the 
feedback due to the propagation of electrons through the lat- 
tice. Metzner and Vollhardt^^ and Miiller-Hartmann^^ noticed 
that the local approximation becomes exact in the limit of in- 
finite spatial dimensions {d °°). The effective site can be 
viewed as correlated atomic problem within a time-dependent 
external field or an effective SIAM.30.3i,49 j^^^ defines the 
self-consistency condition of the DMFT which has been sub- 
ject of two reviews. Within weak-coupling U -perturbation 
theory it could be shown^^ that a local, ^-independent self- 
energy is a good approximation for realistic dimension d ^3 
as corrections due to intersite contributions to the self-energy 
are negligibly small. However, phase transitions remain mean 
field like in DMFT (Refs. 9, 21, and 56) since ^-dependent 
fluctuations are not included in a local approximation. 

The following exact relations for the conduction electron 
Green function Ga{k,z) and the /-electron Green function 
Fa{k,z) can be obtained for the PAM (1): 



Ga{k,z) 



Fa{k,z) = 



Z-Sfa-I.i(k,z) 



z-efa-J:iik,z)' 



\V[ 



-1 



(2a) 



(2b) 



where z is any complex energy off the real axis. Within a 
local approximation such as the DMFT, the ^-dependent /- 
electron self-energy E^(fc,z) is replaced by a ^-independent 



rij(z). From Eq. (2a), one defines a self-energy of the con 



duction electrons via 



Z-£fa-^a{z) ' 



(3) 



which can include a simple k dependence through the hy- 
bridization matrix elements jVp, here taken as constant. For 
such a local self-energy ^^(z), the site-diagonal conduction- 
electron Green function Ga can be written as a Hilbert trans- 
formation 



Ga{z) = ^Y.GaCk,z)^D{z-l.'^,{z)). 
k 

defined for arbitrary complex argument z as 

,Po(e) 



D{z) 



de- 



e 



(4) 



(5) 



where po(fi') is the density of states of the noninteracting con- 
duction electrons. 

The DMFT self-consistency condition states that the site- 
diagonal matrix element of the /-electron Green function of 
the PAM must be equal to Fioctj (z) of an effective site problem 



Mz)^^LFa{lz)^Fioc,aiz) 
k 



(6) 



F\oc,a{z) 



z-efa-^o{z)-^l{z 



with the same local /-electron self-energy {z) for the lat- 
tice and the effective site. This defines the self-consistency 
condition for the functions and A(j(z). 

Given the Green functions Fa{z) and Ga{z), their spectral 
functions determine the local occupation numbers 



(7) 
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^/ d(0 f{(0-ii)lmFa{o}-iQi+), (8a) 



nc = -Y,J dcof{co-n)lmGa{co-iO+), (8b) 



where /(ft)) denotes the Fermi function. Then, the total fill- 
ing per site is given by «tot = «/ + «c- As a matter of conve- 
nience, we will perform an integral transformation such that 
jj. is absorbed into Ef and the band center ; all energies will 
be measured with respect to jj.. For a given lattice filling «tot, 
we have to adjust ji in addition to fulfill Eq. (6). 

Before we will discuss the solution of the effective site, let 
us briefly comment on the implications of the analytical form 
of the conduction-electron self-energy (3). For a Fermi liq- 
uid, the imaginary part of T.-^ vanishes quadratically close to 
the chemical potential for T —> — i.e., lm'L{,{co — /0+) °^ CO^. 
At particle-hole symmetry — i.e., «tot = 2 for a symmetric 
Po(fi)) — E'^(z) diverges like 1 /z leading to an insulator. Away 
from particle-hole symmetry, the denominator remains finite 
and the imaginary part also must have Fermi -liquid properties 
ImE^(a) — /0+) °^ ImE^(a) — /0+) °^ CO^. The real part is very 
large, and therefore, the spectral function of Ga (z) as well as 
Fa{z) samples the high-energy band edges of Pq{co) yielding 
a hybridization gap. This analytic properties must be fulfilled 
by any approximate solution of the DMFT self-consistency 
condition (6). 

The low-temperature physics is determined by the temper- 
ature scale Tiow which is only defined up to a constant factor. 
We will use the renormalization of the Anderson width by the 
quasiparticle spectral weight 



dco 



(9) 



as our choice of such a low-temperature scale Tq oc 7]^^ for 
our numerical analysis. It is related to the mass enhancement 
m* /m = Fo/Tb. We discuss other definitions of such a scale as 
well as the possibility of several low-temperature scales later 
in Sec. VI. 



C. Numerical renormalization group 

The Green function Fioc.a, Eq. (7), of the effective impurity 
problem can be viewed as the /-Green function of an effective 

SIAM,20'21.30.31,50 
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a a •' 

+ Y^j deV^p,sie){dija + fUaE), (10) 



with an energy-dependent hybridization function r(e) = 
^V^Peff(e) = ImA(e — ;0+) describing the coupUng of the / 
electron to a fictitious bath of "conduction electrons" created 
by d^g with density of states (DOS) Peff (e)- The hybridization 
strength V is chosen to be constant, defined via 



TtV" = / der{e), 



and equal to V in the original model (1). 

We accurately solve the Hamiltonian (10) using Wilson's 
NRG.^''-^^ The key ingredient in the NRG is a logarith- 
mic discretization of the continuous bath, controlled by the 
parameter^-' A > 1 . The Hamiltonian is mapped onto a semi- 
infinite chain, where the A^th link represents an exponentially 
decreasing energy scale A^^/^. Using this hierarchy 

of scales the sequence of finite-size Hamiltonians J^n for the 
A^-site chain is solved iteratively, truncating the high-energy 
states at each step to maintain a manageable number of states. 
The reduced basis set of J^^n thus obtained is expected to 
faithfully describe the spectrum of the full Hamiltonian on the 
scale of Dn, corresponding^^ to a temperature Tn ^ Dn from 
which all thermodynamic expectation values are calculated. 
The energy-dependent hybridization function A(z) determines 
the coefficients of the semi-infinite chain. 

The NRG is used to calculate the spectral function of the 
Green function at finite temperatures. We used a slight mod- 
ification of the algorithm^^ for finite temperature Green func- 
tions by Bulla et al.^ with the broadening function^' 

8{co-E) ^ e-'''/4e-[>°g('»/^)/^lV(%/^fe|£|), 

where we choose b ~ 0.6. As usual, the raw NRG Green func- 
tions determine the self-energy {z) by the exact ratio 



pNRG 



iz) 



(12) 



derived via equation of motion technique,*'' where M^^'^ (z) = 
{{faflaf-o\fi))iz). 



D. Modified perturbation tiieory 

The MPT is an approximation for calculating the self- 
energy of the SIAM starting from the following ansatz: 



T.'a{z) = Un'_„ + jT 



y-pa^a\z) 



(13) 



It is based on second-order perturbation theory^^ (SOPT) rel- 
ative to Hartree-Fock solution with U as expansion parameter. 



(2) 

Here, is the second-order contribution to the self-energy 
(choosing Pa =0 reproduces the SOPT), 

j,(2) ^ ^2 fff Pa^(fiJl)p"g(fi>2)p"g(C03) 
JJJ Z-C0i + 0>2-0>i 

x{/(«i)[l-/(«2)]/(«3) 
+ [1 -/(t0i)]/(C02)[l - f{(03)]} d(Oi da)2d(03, (14) 

with the Hartree-Fock /-electron spectral functions of the ef- 
fective SIAM: 



(11) p^P(a)) = --Im 



1 



CO+iQ+ +jl-efa- UnL„ - ^a{(O+i0+) 

(15) 

The parameter j3(j is constructed such that the exactly solvable 
atomic limit V = as well as the first four spectral moments 
are correctly reproduced. ^^-^^ The Hartree-Fock / occupation 
must equal the full / occupation determining the effective 
Hartree-Fock chemical potential ji. 



III. RESULTS FOR THE SINGLE-PARTICLE DYNAMICS 

All transport calculations rely on the results for the single- 
particle Green functions of the PAM. We use a Gaussian 
model density of states for the unperturbed conduction- 
electron system — i.e., po(fi)) = exp[-(co/f*)^/2]/(f*\/27r) — 
which is appropriate for a li-dimensional hypercubic lattice 
in the limit d ^ o°}'^ In the following, we measure energies 
in units of To ~ TtV^pQ{Q) — i.e., the hybridization is fixed at 
= 2fTo/%/27r — and choose \/2f* = IOFq. As impurity 
solvers, we use the NRG as well as the MPT. The chemical 
potential /i has to be determined self-consistently for a given 
total number of electrons, «tot, per site; in the following fig- 
ures, energies (frequencies) are measured relative to the chem- 
ical potential pi — i.e., we have jU = and the band center is 
shifted accordingly. 



A. / -electron spectral function and self -energy at finite 
temperature T 

Using the NRG the calculation of the spectral functions 
for finite T is numerically very challenging. In contrast to 
the spectral functions^^'^' at T = 0, where only the ground- 
state excitations are relevant, all excitations contribute ac- 
cording to their statistical weight, but the NRG only pro- 
vides spectral information up to an energy scale con °= A^^l^. 
While the number of NRG iterations, A^, is in principal arbi- 
trary and arbitrary small excitations could be resolved, Wil- 
son has already pointed out^^ that the lowest energy scale cOn 
should be identified with the temperature T for which ther- 
modynamical expectation values are calculated. It is obvi- 
ous from the Lehmann representation of the spectral func- 
tion that eigenstates of the Hamiltonian with eigenenergies 
Em < T contribute equally to the spectral functions. There- 
fore, we stop the NRG iteration when ~ T and must inter- 
polate for frequencies \(o\ < (On- The technical details of our 
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FIG. 1: /-DOS P/((b) calculated with DMFT-NRG for U /Tq = 5, 
Ef — ec = — 2.5ro, chemical potential ^ = 0, and a filling ;7tot =1.6 
at several finite temperatures T /Tq. NRG parameters: number of 
retained NRG states, N, = 1500, A = 1.6, S/Fq = 10"^ 




FIG. 2: /-DOS p/(co) calculated with DMFT-NRG for U /Tq = 10, 
£f — £c = — SFo, chemical potential 11 = 0, filling «tot = 1.6, and 
finite temperature T = O.OOOBFq. A three-peak structure with peaks 
at £f, and £f + U is visible. In the inset the transition from a 
hybridization gap for U = (exact result) to a pseudogap for U /Fq e 
{3,6,50} is shown. NRG parameters: as in Fig. 1. 



finite-temperature algorithm are described in the Appendix of 
Ref. 59; other possible approaches are found in Refs. 60 and 
63,64,65. 

An exemplary series of /-spectral functions P/(ft)) calcu- 
lated with DMFT-NRG for different temperatures T and fixed 
U = STq is shown in Fig. 1 for the metallic regime with a 
filling «tot = 1 .6 and Ef — Ec = — 2.5ro. A hybridization pseu- 
dogap develops close to the chemical potential {1 = for T 
decreasing below a characteristic temperature Tjow. The mod- 
erate value of U leads to a spectrum where the high-energy 
charge excitations are not well separated from the resonance 
close to the chemical potential. 

The typical structure of the /-spectral function Pf{(o) cal- 
culated with DMFT-NRG for a larger value of U /Fq = 10, 
Ef — Ec = — SFo, at a temperature T ~ O.OOOSFo well below 
Tjow is depicted in Fig. 2 (a detailed discussion of possible 
estimates for Tjow is found in Sec. VI). A pronounced peak 
structure with a pseudogap dominates the low-energy part of 
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FIG. 3: Comparison between the /-DOS P/(tt)) calculated with 
DMFT-NRG and DMFT-MPT for (//Fq G { 10, 5,3,1} and £f~Ec = 
—U /2 at finite temperatures T /Tq as indicated. Here, the chemical 
potential is /i = with filling «tot = 1-6. NRG parameters: as in 
Fig. 1. 



the spectrum in the vicinity of the chemical potential similar to 
the one shown in Fig. 1. In addition, we observe two shallow 
high-energy peaks, one at Ef below /i, one at + [/ which 
corresponds to double occupancy of the / levels. It is always 
correctly positioned by the NRG, independent of the value of 
U , but with a linewidth too large due to the NRG broadening 
procedure (see Refs. 59 and 61 for details). 

In the inset of Fig. 2 the development of the hybridization 
gap of the /-spectral function P/(o) with increasing U > Q 
is shown and compared with the exact result for U = Q. For 
U = 0, we obtain from an uncorrelated hybridized band a hy- 
bridization gap above = of a width oc V . It arises from 
the divergence of the real part of {z^ Ef) given in Eq. (3): 
the high-energy part of the free density of states po defines the 
shape of the gap via the exact relation between and Fa, 



Faiz) 



y2 



{l+L%iz)D[z~L,{z)]}. (16) 



For a DOS Po{(o) of Gaussian shape the hybridization gap is, 
strictly speaking, only a pseudogap even atU = 0. A real gap 
only arises for densities of states with well-defined sharp band 
edges such as a semielliptical or a true finite-c/ DOS. However, 
with increasing U the gap will always evolve into a pseudo- 
gap anyway due to the finite lifetime of the quasiparticles, re- 
flected in the growing imaginary part of E^(a)). This effect is 
clearly seen in the inset to Fig. 2. The width of the pseudogap 
and of the whole Kondo resonance peak structure is decreas- 
ing with increasing U stemming from the the reduction of the 
quasiparticle weight at the chemical potential. 

A comparison of the /-spectral functions for finite T ob- 
tained in DMFT-NRG and DMFT-MPT for Ef ~ Ec ^ -U /2 
and «tot = 1.6 (metallic case) and four different values of U 
is displayed in Fig. 3. Again, only a pseudogap is found at 
finite U which is narrowing with increasing U. As expected, 
the DMFT-MPT curves agree very well with the DMFT-NRG 
graphs for small U. This is not surprising, because the MPT is 
based on U -perturbation theory and must, therefore, become 
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FIG. 4: (Color online) Comparison between the self-energy {(o) 
in DMFT-NRG and DMFT-MPT for U /Tq e {3,5, 10}. The in- 
set shows a close-up of the vicinity of = 0. All parame- 
ters are as in Fig. 3, in particular finite temperatures T /Tq e 
{0.03,0.019,0.0029}. 
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FIG. 5: Influence of the shift 8 on lmZ''(^i,r) for (//Fo = 10, £/- 
Ec = -U/2, and «tot = 1.6. The ratio r = y2/[e^ + ReZ/(0)]2 (right 
inset) is almost constant for low T, here r ~ 14 1. This leads to 
the enhancement of ImS'^'(/x, T) for low T compared to ImE^ (Mi ^) 
as shown in the left inset. NRG parameters as in Fig. 1 . 



correct for sufficiently small U . Already for U /Yq = 5, how- 
ever, one observes deviations between the MPT and NRG ap- 
proaches. For U /To = 10 {U j^/lt* = 1) the width of the reso- 
nance peak at the Fermi energy is obviously too large in MPT; 
i.e., quantitatively there is a strong overestimation of the low- 
energy scale within DMFT-MPT. For large U the MPT cannot 
reproduce the correct Kondo temperature scale, as the physics 
in this regime is driven by the spin-flip scattering J propor- 
tional to 1/f/ and is, therefore, nonperturbative in JJ . 

Figure 4 shows the frequency dependence of the (imagi- 
nary part of the) /-electron self-energy obtained in DMFT- 
NRG and DMFT-MPT, corresponding to Fig. 3. Here, the 
imaginary part of the self-energy is finite even at o = 0, and 
away from the chemical potential one observes a quadratic 
(o) — yif' behavior as expected for Fermi liquids. The finite 
value at o = /i = is not only due to the finite T , but it has 
also a contribution from a finite imaginary part 5 (i.e., an ad- 
ditional Lorentzian broadening), which we had to introduce 
for numerical reasons. As we observed numerical instabilities 
in the self-consistency equations (6) and (7) as well as inac- 
curacies in the local self-energy (12) in NRG, we solved the 
set of equations (6) and (7) in the complex plain at a finite 
shift 5 /Fo = 10^^ away from the real axis to obtain stable nu- 
merical solutions. Physically this artificial broadening 5 can 
be interpreted to simulate the effects of impurity scattering, 
yielding a finite lifetime corresponding to a finite imaginary 
part of an additional "disorder" self-energy. But one also sees 
from Fig. 4 that — in spite of the same small imaginary part 
5 — the MPT result for ImE^(O) is systematically smaller than 
the NRG result, even for small V . This is an indication that the 
finite value in NRG is not only caused by the temperature off- 
set and 5 but also by the additional Gaussian broadening^^'^' 
described in Sec. II C and additional numerical errors stem- 
ming from the relation (12) to determine the self-energy in 
NRG at finite temperature. 



B. Temperature dependence of band-electron self-energy 

Since the conduction-electron self-energy E'^(a),r) 
strongly influences the transport properties of the PAM, we 
study its imaginary part at to = — i.e., the conduction- 
electron scattering rate. According to Eq. (3), it is given 
by 



ImE',(0) = 



|VpImE{(0) 



[e/a+ReE'^(0)]2 + [ImEi(0)]2 



(17) 



Analytically, we can distinguish two cases. 

(i) If |e/a + Rela(0)| < |ImE^(0)|, the conduction- 
electron scattering rate is reciprocal proportional to the /- 
electron scattering rate — i.e., ImE^(O) w y^/ImE^(0). This 
is the case for the Kondo insulator regime, where Zf^ + 

ReE^(O) ^0. 

(ii) Only when |e/a + ReE{(0)| > |ImE{(0)| is the de- 
nominator of Eq. (17) dominated by the real part [e^cr + 
ReE'^(O)]^ and the c-scattering rate becomes proportional to 
|ImE^(0)|. 

We show the temperature dependence of the c-electron 
scattering rate and its dependence on 5 for a fixed value 
of U /Fo = 10 in Fig. 5. While the /-electron scattering 
rate remains very small for low T (cf. left inset of Fig. 5), 
the c-electron scattering rate has a much larger finite T 
value. This originates from the small but finite imagi- 
nary part of the conduction-electron self-energy (17) as a 
consequence of the small (artificial) imaginary part 5, here 
S/Fo e {10^2 5 X 10"^ 10"^ 10"'^} and the numerical error 
in the ratio of two Hilbert transformed spectral functions in 
Eq. (12). The small finite |ImE{-(0)| is enhanced if the ratio 
r /[Efa +ReE'^(0)]^ is larger than 1, as it is usually the 
case since the hybridization gap is very close to the chemical 
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FIG. 6: (Color online) Contour plot of the total energy-dependent 
density of state p(e,co) = p/(e, cu) + Pc.(e, to) for U /Tq = S, £f~ 
Ec = -U/2, «tot = 1.6, /i = 0, T/Fq = 3 X 10""*. NRG parameters: 
as in Fig. 1. 

potential, and —£fa — ReE^(O) is an estimate for its location. 
As one can see from the right inset of Fig. 5, we indeed have 
a ratio r > 14 > 1. Therefore, we can estimate the c-electron 
scattering rate by 

lmi:^(0) - 14ImE{(0) 

for the chosen parameters U /To = 10 and «tot = 1.6, which 
explains why the artificial finite imaginary part is even 
more important and pronounced for ImE'o^(O) than it is for 
Iml'^(O). 

C. Renormalized band structure 

In the previous section. Sec. Ill A, we presented the local, 
or ^-summed spectral function. In the Fermi-liquid phase of 
the model, new quasiparticles are formed as a mixture of / 
and c degrees of freedom. The resulting renormalized band 
structure can be visualized by plotting the energy £ = and 
frequency dependent density of states p(e,a)) = Pf{£,co) + 
Pc{£,(o), where 

Pf{e,(o) = lmFcr{e,(o-iQ^)/n, (18a) 
Pci£,co)=lmGa{e,co-iO+)/7t, (18b) 

as a two-dimensional (2D) color (online) contour plot dis- 
played in Fig. 6. For «tot = 1-6, the chemical potential jj. =0 
lies at the top of the lower hybridized band. For V = 0, 
we would have a sharp line at o = e — with purely c 
character, while two lines with fractional weight of 1/2 at 
CO = Sf = Ec - U /2 and co = Ef + U = Ec + U/2 with / char- 
acter could be found. At finite V, the / electrons become part 
of the Fermi volume. Indeed, the almost dispersionless and 
therefore heavy quasiparticles close to the chemical potential 
have mainly / character. The second band is located above 
the chemical potential. We note — not shown here — that this 
renormalized band picture remains valid even well above the 
low-temperature scale. The dark colored peaks of p (e, co) are 



broadened by the increasing imaginary part of the self-energy 
but can still be traced by the zeros of the real part of the re- 
ciprocal Green functions Ga{E,Co) and Fa{£,Co). This ob- 
servation turns out the be very important to understand the 
midinfrared peak in the optical conductivity. ^^'^^ 



IV. TRANSPORT THEORY 

To describe the electronic transport within the PAM we start 
from the standard relations for the generalized transport coef- 
ficients, according to which the electrical current density J and 
the heat current density q depend linearly on the electric field 
E and the temperature gradient VT: 

J = LnE + Ln{-\^yT), (19a) 
q^L2iE+L22{-^VT). (19b) 

All coefficients are calculated within the linear response ap- 
proach, starting from similar Kubo formulas. ^^-^^ For symme- 
try reasons, L12 = Lai holds. 

For example, the real part of the frequency dependent (op- 
tical) conductivity tensor^^-^^'™ g{co) = Lii(a)) is related to 
the current-current correlation function and written as 

<ya^((0) = -^Im(0„|jJ))(a) + /0+), (20) 

where Vo = is the volume of the unit cell and counts 
the number of lattice sites. It has been shown^' that the cur- 
rent operator of the PAM has two contributions; a conduction- 
electron part and a part proportional to W^. The / electrons 
do not appear in the current, since they do not disperse. For 
a ^-independent hybridization, only the conduction electrons 
carry the electrical and heat currents: 

/^^E^L^Ja' (21) 
ka 

where = ^V^gj is the group velocity. Hence, the current- 
susceptibility tensor {{j\j^)){z) is connected to the particle- 
hole Green function 

imm^e' L %-1«l-.j4.,^r'.'))W- (22) 
aa'kk' 

In a cubic crystal, the conductivity is isotropic: Gap (co) = 
G{co)i. From now on, we will consider only the xx- 
componentof the conductivity g{co) = o;,x(fi)). 

In general, the full two-particle Green function 
((<^7 <^T,TkT involves vertex corrections which 

reflect residual particle-particle interactions.^^ However, in 
the limit d ^ oa \i was shown that current operator vertex 
corrections vanish. ^^-^^ Thus, it is consistent with the DMFT 
assumption of a ^-independent self-energy that these vertex 
corrections vanish for any lattice model of correlated electron 
systems. For the special case of a local approximation for the 
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PAM this was already shown in Refs. 41, 42, and 74, as for 
symmetry reasons 

£vjH^,|2^j(z + a))Gj(z)=0. 

k 



Therefore, we obtain 



((y,|/:))(a) + /0+) = f ^ ) / ./a)7(«')P.(ej,«') 



X [Ga(^,a)' + ffl + /0+) 

+ Ga(^,«'-fi)-/0+)]. 



(23) 



Within the DMFT, the lattice one-particle Green function de- 
pends only on the (complex) energy z and bare band disper- 
sion £j: Ga{k,z) = Gc7(ej,z). Then 

]5^L(^) ^^^l)^ |_^^/epo(e)A(£), (24) 



with 



(25) 



The thermoelectric power S is defined as the proportional- 
ity constant between an applied temperature gradient and the 
measured voltage drop in the absence of a current flow. The 
Peltier coefficient given by the ratio of heat and electrical cur- 
rent is related to the thermoelectric power by n = TS. Note 
that \e\L\2lkBTL\\ is dimensionless and kB/\e\ « 86 /iV/K. 
Therefore, the thermoelectric power is given in absolute units; 
only the scale of temperature axis must be fixed by experi- 
ment. 

If we assume one electron per unit cell of the volume a^. 



10 m, the resistivity p ^ a has the natural unit 
C7q"' = ha/e^ w 41 jUi2 cm. 



(31) 



Similarly, if we assume Tq = 100 meV, the thermal conduc- 
tivity is given in units of 



ha 



' Q.2\ W/(Kcm) 



(32) 



Finally, the dimensionless figure of merit is defined as 



(33) 



which measures the efficiency of a thermoelectric material. 



Po(e) has been evaluated approximately in large dimensions'^ 

as 



Po(e) = ^Po(e) + 0(l/«f2) 



(26) 



on a hyper-cubic lattice. Then, Eq. (23) can be reduced to a 
sum of Hilbert transforms, defined in Eq. (5). 

By taking the limit o — > in Eq. (23), the static conductiv- 
ity a ~ L\i is obtained, and we get for the generalized trans- 
port coefficients: 



Ln^j^j_yf'{(o)]x[(o)d(o, (27a) 
Ln = ^J_J-f'{(o)]{co-^)Tico)d(0, (27b) 

L2j = ^[ [-f'{(o)]i(0-^fT{(o)d(o. (27c) 

na j-oo 

Here t{(o) represents a generalized relaxation time defined as 
T{(o)^—{t*f po{e)p^ie,co)de (28) 



and /'(ft)) is the derivative of the Fermi function. Only in the 
Fermi-liquid regime does one have T(ftj) °^ l/ImEc(ftJ + /0+) 
and the linearized Boltzmann transport theory is recovered. 

The thermal conductivity k and the thermoelectric power S 
are then given by^^'^" 



TLn \e\kBTLi 



1 
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(29) 
(30) 



V. RESULTS FOR THE TRANSPORT PROPERTIES 
A. Resistivity 

The (static or direct current) resistivity p{T) is obtained 
from the reciprocal of (j(r), Eq. (27a). To achieve bet- 
ter convergence of the DMFT equation at low temperatures 
again a finite offset 5 > from the real axis was intro- 
duced. As discussed in Sec. Ill A, the NRG provides ac- 
curate spectral functions only for frequencies co > T. In 
the linear-response transport-integrals (27), however, only 
the low-frequency spectral information in the Fermi window 
|ftj| < 2T contributes. This intrinsic property of the NRG 
makes it very difficult to obtain reliable spectral information 
for low frequencies at finite temperatures needed for calculat- 
ing transport properties. 

In Fig. 7, the resistivity is displayed for fixed filling «tot = 
1.6, U/Fo = 10, Vlt* = lOFo and various values of the shift 
5 as in Fig. 5. Physically, a finite 8 can be interpreted as 
simulating the effects of lattice defects and/or of a small but 
finite concentration c of nonmagnetic impurities giving rise to 
a finite "impurity self-energy" imaginary part (inverse scatter- 
ing time) S ^ c, which leads to a finite residual resistivity for 
r ^ 0. Figure 7 shows that the residual resistivity p{T — Q) is, 
in fact, increasing with increasing 5. As discussed already in 
Sec. IIIB, the finite imaginary part of the conduction-electron 
self-energy ( 17) is much enhanced compared to the imaginary 
part of the /-electron self-energy, which explains the strong 
dependence on 5 seen in Fig. 7. However, even in the limit 
5 a finite p (r = 0) is obtained for T ^ in DMFT-NRG. 
Not only the finite 8 but also the additional Gaussian broaden- 
ing and the limited accuracy of the NRG in the regime a <T 
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FIG. 7: Influence of the shift 5 on the static resistivity in DMFT- 
NRG for U /To = 10, £/ — Ec = — f//2, and ;itot = 1 -6. All parameters 
are as in the corresponding Fig. 5. 
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FIG. 8: Resistivity p{T) obtained in DMFT-NRG for different num- 
ber of kept states A'^ and fixed values of S/Tq = 10^^ and A = 1.6. 
Other parameters as in Fig. 7. 



causes the finite self-energy imaginary part and the finite p (0). 
We reproduce the typical metallic HFS behavior'^'^^-'^ within 
our DMFT-NRG treatment: a resistivity increasing with in- 
creasing T for low T, a maximum of the order 100 jLiii cm at a 
characteristic temperature Tmax (about 10-200 K), and a p{T) 
(logarithmically) decreasing with increasing T for T > Tmax- 

For the same parameters as in Fig. 7 the dependence of 
p{T) on the NRG parameter Ns (number of states kept) is 
shown in Fig. 8. One observes that for Ns > 1500, p{T) no 
longer depends on Ns- Therefore, all further DMFT-NRG cal- 
culations have been performed for this Ns = 1500. 

For fixed S/Tq ~ 10^^, a total filling ntot = 1.6, a differ- 
ence between the band center and Ef of Ef — Ec = —U /2 
and several small values of U the DMFT-NRG and DMFT- 
MPT results for p{T) are compared in Fig. 9. One observes 
that in the high-temperature regime both results agree and that 
the overall agreement becomes the better the smaller U is, as 
expected already from our discussion in Sec. Ill A, Fig. 3. But 
for the same small 5 the residual resistivity p(0) in DMFT- 
MPT is much smaller than in DMFT-NRG. In DMFT-MPT 



FIG. 9: (Color online) Comparison between the resistivity p{T) in 
DMFT-NRG (solid lines) and DMFT-MPT (dashed lines), for dif- 
ferent U , £f — Ec = —U /2, and fixed total filling «tot = 1.6 (metallic 
case). NRG Parameters: as in Fig. 1. 



the finite p (0) is solely determined by the finite 5 whereas in 
DMFT-NRG the additional broadening and numerical inaccu- 
racies at very small (O contribute. 

For the same parameter set, but higher values of the corre- 
lation U, DMFT-NRG resistivity results are plotted in Fig. 10. 
The temperature Tmax, at which the resistivity has its maxi- 
mum, is shifted to lower values for increasing U; as discussed 
in Sec. VI one obtains an exponential dependency of TJnax 
on U. For sufficiently strong U the peak height at Tmax is 
nearly U independent. This is seen once more from Fig. 10(b), 
which shows the scaling properties plotting the resistivity 
p (r) — p {To /lOO) versus T /Tq [for the low-temperature scale 
To see Eq. (9) and Sec. VI]. While for U/Fq < 7, where high- 
and low-temperature scales are not very well separated, the 
maximum of the resistivity and the peak height show a U de- 
pendence, we reach a universality regime for large U. Note 
that To is of the order of the position Tmax of the maximum of 
the resistivity while the T^ behavior of the resistivity is only 
observed below the coherent scale Tcor which is two orders of 
magnitude smaller 

The dependence of the resistivity p{T) on the total oc- 
cupation «tot is shown in Fig. 11 for fixed U /To = 10 and 
Ef — Ec: = —U /2. We see that the (artificial) residual resis- 
tivity p (0) increases whereas the maximum temperature Zmax 
decreases when approaching the Kondo insulator regime at 
«tot = 2 (symmetric PAM). In this case p{T) saturates for 
T < 10^^, which is typical for Kondo insulators.^ Experimen- 
tally, this is due to impurities; in our calculations it is a conse- 
quence of the finite-energy imaginary part 8. For fixed «tot = 2 
and different values of Ef — Ec, p{T) is shown in Fig. 12; obvi- 
ously, the Kondo insulator behavior is not only obtained in the 
symmetric case but always for «tot = 2 in accordance with the 
Luttinger theorem. The (artificial) finite p (0) depends hardly 
on £/ — Ec but only on the finite 8. 
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FIG. 10: (Color online) Resistivity p{T) as function of T for dif- 
ferent U, Ef — Ec = —U /2, and fixed total filling ntot = 1.6. The 
calculations were done in DMFT-NRG, thus continuing the NRG se- 
ries of Fig. 9 to higher U. NRG Parameters: as in Fig. 1. In (b) the 
same data is plotted as p{T) — p (To/ 100) +5 fl^l cm versus T/Tq 
with Tq as in Eq. (9). 




FIG. 12: Resistivity p{T) for different £/-£(-, S/Tq = 10"^ 
U /To = 10, ;itot = 2 (Kondo insulators). NRG parameters: as in 
Fig. 1. 
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FIG. 13: Optical conductivity for four different temperatures T : 
0.01,0.1,1, low, "tot = 1-6, and a fixed U /Fq = 10, e/-ec 
—U /2. NRG parameters: as in Fig. 1. 




FIG. 11: Resistivity p{T) for different total filhngs «tot, U /To = 10, 
£f-Ec = -U /2, 5/ro = 10"\ NRG parameters: as in Fig. 1. 



B. Optical conductivity 

The optical conductivity (j(ft)) obtained in the metallic case 
for «tot = 1 .6 is shown in Fig. 13 for different temperatures T. 
For very low T one observes a Drude peak at low frequen- 
cies and, in addition, a "midinfrared" peak at finite frequency 
CO ^ Vf*To- With increasing T the Drude peak quickly de- 
creases (in accordance with the strong increase of the static 
resistivity), whereas the midinfrared peak remains nearly un- 
changed, and the remainders of the Drude peak and the mid- 
infrared peak merge into a broad structure in the high temper- 
ature regime. 

Figure 14 shows CJ(ft)) for a Kondo insulator, where for 
very low T the Drude peak is absent and only the midinfrared 
peak is present due to interband transitions. With increas- 
ing T, however, the gap disappears, and the metallic heavy- 
fermion behavior at high temperatures is recovered. 

This behavior can easily be understood in terms of the 
renormalized bands discussed in Sec. Ill and shown in Fig. 6. 
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FIG. 14: Optical conductivity for three different temperatures T, 
(itot = 2 (i.e., Kondo insulators), and a fixed U /Fq = 10, £f — £c = 
—U /2. NRG parameters: as in Fig. 1. 



With the assumption of a constant density of states Po(m), 
cj(a)) can be approximated by 



a{co) _ ipoip) 



2co 



dco'[fico')-f{co' + co)] 



(0 + T.'{co' - i5) - T.'{co' + CO + i5)' 

(34) 

The optical conductivity (J (ft)) has a maximum when the real 
part of the denominator on the right-hand side of Eq. (34) 
vanishes. This is obviously the case for o ^ for metal- 
lic HFSs in the Fermi-liquid regime, which yields the Drude 
peak. An additional maximum is found in the vicinity of the 
minimal intraband transition energy close to the chemical po- 
tential, which yields the midinfrared peak. An estimate for its 
position is the smallest interband transition energy for ^ = 
given by Aopt °^ Vt*Jo which fits perfectly with the observed 
midinfrared peak of a large variety of different heavy-fermion 
compounds.^" Its temperature dependence is weak since the 
spectral weight moved from lower to higher excitation ener- 
gies is small. On the other hand, the midinfrared peak width 
depicted in Fig. 13 and reported in Ref. 76 comes out much 
narrower than the experimentally observed peaks. This is 
expected because in real materials, all Hund's-rule multiplets 
contribute to the scattering of the conduction electrons in the 
midinfrared region. Our simplified model, however, only con- 
tains the lowest Hund's-rule doublet. The total optical re- 
sponse is furnished by a superposition of all scattering res- 
onances which leads to an effectively broad resonance as ob- 
served experimentally. 



C. Thermoelectricity 

The thermoelectric power S{T) measures the ratio between 
electrical and heat cuiTent divided by the temperature, and its 
sign is related to the integrated particle-hole asymmetry rela- 
tive to the chemical potential. In Fig. 15, S{T) is plotted for 
Ef — = —U /2 and fixed total filling for a variety of differ- 
ent values of U [for the low-temperature scale Tq see Eq. (9) 
and Sec. VI]. We obtain very large absolute values for S{T), 
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FIG. 15: (Color online) Thermoelectric power for various values of 
U and fixed total occupation ntot = 1.6 and Ej — £c = —U /2. Panel 
(b) shows the same date but on a rescaled axis T /Tq with Tq as in 
Eq. (9). All Parameters: as in Fig. 10. 



of the magnitude 50 juV/K up to 150 /iV/K; note that the 
thermoelectric power is obtained in absolute units, as already 
mentioned in Sec. IV. Similar to the resistivity, the thermo- 
electric power exhibits a low-temperature peak which is corre- 
lated with the maximum of the resistivity which is an analyti- 
cal consequence of Eq. (29). In addition, we observe a second 
extremum moving to higher temperatures with increasing U 
which results from the charge fluctuations on the energy scale 

Ef-ll. 

It is very interesting to note that for U /Yq > 8 the ther- 
moelectric power changes sign in an intermediate-temperature 
regime. Several extrema and similar sign changes were found 
experimentally for several HFS materials, in particular below 
80 K for CeCu2.2Si2."'^^ As we have discussed in Sec. Ill, 
the details of the spectrum close to the hybridization gap de- 
pend on the free single-particle DOS po(co). In addition, our 
model does not include other Hund's-rule multiplets. There- 
fore, we only can make a qualitative comparison with exper- 
iments. Remarkably, however, we find a sign change of the 
thermoelectric power for large U and a large absolute value of 
5(7") comparable to the experimental results. 

The thermoelectric power 5(7") for different fillings «tot, 
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FIG. 16: Thermoelectric power for a fixed U /Fq = 10 and various 
occupations ntot = 1-4, 1.6, 1.9, 1.95,2 and £f — £c = —U /2. All pa- 
rameters: as in Fig. 11. 




FIG. 17: Thermoelectric power for U /Fq =10, /itot = 2 (Kondo insu- 
lator regime) and different values for £f (or ey — £c). All parameters: 
as in Fig. 12. 



U /To = 10, and e/ — = —U/2 is depicted in Fig. 16. For 
«tot = 2 the thermoelectric power vanishes for all tempera- 
tures because of particle-hole symmetry. Otherwise S{T) ex- 
hibits a peak with increasing height when approaching the 
Kondo insulator regime ntot ^ 2. For the Kondo insula- 
tor situation ntot = 2 and various £/ — the thermoelectric 
power S{T) is shown in Fig. 17. Only in the symmetric case 
Ef = —U /2 — — SFo does S{T) vanish due to particle hole 
symmetry; away from the symmetric case, the peaks reach 
extremely large absolute values of S up to 200 /zV/K. This 
shows that Kondo insulators could be candidates for thermo- 
electric applications. 

However, the figure of merit Z = S'^a/fC, or the dimension- 
less value ZT and not only the Seebeck coefficient S(T) is 
of relevance for efficient thermoelectric cooling, where k is 
the thermal conductivity. Since approaching the Kondo insu- 
lator reduces the conductivity a as well as shifting the peak 
to very low temperatures the nominator of ZT will be rather 
small. Fig. 18 shows the electronic part of the thermal con- 
ductivity k{T). As expected k{T) is rapidly increasing with 



FIG. 18: Electronic part of thermal conductivity for U /Fy = 10, Ey 
£(■ = —U /2. All parameters: as in Figs. 1 1 and 16. 
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FIG. 19: Figure of merit for U /Tq = 10, ef-£c = -U/2. All pa- 
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increasing T and the overall behavior is similar as that ob- 
served experimentally. 

The resulting temperature dependence of the figure of merit 
is shown in Fig. 19. According to our calculations a figure of 
merit larger than 1 can be obtained in HFSs in certain low- 
temperature regimes, showing that, in fact, they are candidates 
for thermoelectric applications. 



VI. LOW-TEMPERATURE SCALES 

There has been a long debate whether low-temperature 
properties are governed by one or more low-energy scales in 
the paramagnetic phase of the PAM.^^ '*^'^'''^^'^° In the SIAM, 
the screening of the local moments occurs on the scale of 
Kondo temperature Tk where the system crosses over from the 
local-moment to the strong-coupling fixed point when low- 
ering the temperature. But Nozieres has argued that for the 
lattice with a macroscopic concentration of / moments there 
is an "exhaustion" of the conduction electrons^^'^" available 
for screening of those / moments. However, it is question- 
able whether the concept of individually screened indepen- 
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dent single impurities still holds in a lattice, because one con- 
duction electron could rather contribute to the screening of 
several local moments. In this context, it has been suggested 
that two different low-temperature scales exist in the param- 
agnetic phase of the PAM: namely, the single-impurity Kondo 
scale Tk and a strongly reduced lattice scale Tjow <C Tk at 
which the screening of the local moments happens. It was 
argued that at Tk only «scr = Pc{0)Tk conduction electrons 
are available to screen itf local / electrons, where Pc(0) is 
the conduction-electron DOS at the chemical potential. How- 
ever, in Nozieres's work^^ only phenomenological estimates 
and no precise mathematical definitions for different temper- 
ature scales are given. 

It is tempting to identify the possibly existing additional 
scale 7]ow with the "coherence temperature" — i.e., the scale on 
which the heavy quasiparticles are formed and/or on which the 
behavior in the resistivity is observable. Experimentally, 
however, no evidence was found for such a large separation 
of energy scales in heavy-fermion compounds.''^ On the con- 
trary, substitution of Ce by La in Cei-xLarCug exhibits nearly 
perfect scaling of the Ce contribution to the specific heat'^ 
from the diluted impurity situation (jc w 1) to the lattice case 
{x = 0) indicating that there is only one relevant energy scale. 
Recent theoretical investigations by Vidhyadhiraja and Logan 
gave evidence for such a single low-temperature scale,^^'^^'^^ 
using Fermi-liquid assumptions for the /-electron self-energy. 
Their arguments are consistent with the notion that a new 
thermodynamically relevant energy scale can only occur as 
crossover scale from competing low-temperature fixed points. 
Our NRG calculations for the effective site, moreover, gave 
no evidence for new fixed points of the SIAM other than 
those already discussed in Ref. 57. In fact, we always find a 
strong-coupling fixed point independent of the filling connect- 
ing Kondo insulators adiabatically to metallic heavy-fermion 
systems. 

In the literature, different notions of what is implied by 
the term "low-temperature scale" have led to a controversy 
which could have been resolved by stating clearly the pre- 
cise mathematical definitions of such distinct scales. The 
low-temperature scale Tiow allows to map any temperature- 
dependent property P{T) over a limited low-temperature 
range^"* onto a corresponding dimensionless universal func- 
tion p{x = T /T\ov/) where all details of the initial microscopic 
Hamiltonian enter this single scale Tlow From a renormal- 
ization group perspective, it describes the crossover from one 
fixed point to the T = fixed point. If such a quantity P has 
distinct analytical behavior in different temperature intervals, 
it is tempting to assign additional low-temperature scales to 
identify such different regimes. If, however, a universal func- 
tion p{x) can be found such that P{T) collapses onto this 
curve for all low-temperature regimes, these distinct scales 
contain no additional physical information. The confusion 
in the literature arises by ignoring that these admittedly dif- 
ferent low-temperature scales might ultimately be connected 
and reflect the same physical phenomena. If, however, differ- 
ent scaling regimes can be identified using two or more low- 
temperature scales, it should correspond to additional fixed 
points of the Hamiltonian. 




0.5 1 1.5 2 2.5 3 3.5 
1/g 



FIG. 20: The four different low-temperature scales plotted versus 
the reciprocal dimensionless coupling constant 1/g [cf. Eq. (35)]. 
The inset shows the same data as function of U. The stars mark 
the maximum of the resistivity Tmax/ro the crosses the screening 
of the spin moment of the effective site. The squares indicate Tp = 
1 / \/ooA, while the mass enhancement factor is plotted by the circles. 
Parameters: £f — £c = —U /2, ntot = 1-6, number of retained NRG 
states N, = 1500, A = 1.6, S/Tq = 10"^ 

A. Definition of different low-temperature scales 

From our data, we have defined /oMr different possible low- 
temperature scales characterizing different physical phenom- 
ena in order to identify more than one low-temperature scale. 
These scales are (i) the low-temperature scale^^ Tq = Fom/m*, 
Eq. (9), defined by the quasiparticle spectral weight, measur- 
ing the reciprocal mass enhancement m/m* for T ^0, 



(ii) the A coefficient of the resistivity p{T) = po +AT^ valid 
in the coherent Fermi-liquid regime, Tp ~ (ctoA)^'/^, where 
' is the natural unit of the resistivity given in Eq. (31); (iii) 
the position of the maximum of the resistivity curve T^sf. for 
metallic heavy-fermion systems; and (iv) the screening tem- 
perature T* of the impurity moment^^ (5?)imp(r*) = 0.05. 

Other possible scales such as the width of the hybridiza- 
tion gap are not very well defined since the position and depth 
of the gap as well as the shape of the spectral function are 
strongly dependent on the Fermi volume. In the following, we 
will show that all scales are connected by a fixed ratio when 
varying the interaction strength JJ for fixed occupancy ntot or 
Fermi volume. 

The mass enhancement factor m* /m as well as the A coef- 
ficient of the resistivity contains information on the low-lying 
excitations of the lattice Fermi-liquid fixed point as pointed 
out by Logan and collaborators.^^ Both should be proportional 
to a low-temperature scale Tiow We contrast this with two es- 
timates of a possible second low-temperature scale. It was 
argued that this scale governs the crossover from a regime of 
incoherent scattering of the conduction electrons to a coher- 
ent Fermi-liquid ground state of the lattice. This truly new 
scale would mark the deviation of a periodic Anderson model 
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FIG. 21: (Color online) The temperature dependence of the effective 
moment (5?)in)p(r) versus T /Tq for various values of U and a fixed 
filling of ntot = 1-6. Parameters: as in Fig. 20. 



FIG. 22: (Color online) The temperature-dependent Kondo tempera- 
ture Tf[{T) of the effective site for various U plots versus the scaling 
variable x=T /Tq for fixed filling. Parameters : as in Fig. 2 1 . 



from the properties of the SIAM. We take the position of the 
maximum of the resistivity Tmax as one clear indicator of such 
a scale as well as the screening of the impurity effective mo- 
ments which defines the temperature T* . However, we like to 
point out that our definitions of the low-temperature scales do 
not imply an energy hierarchy Tq < Tmax- 

B. Discussion of the scales 

Using Wilson's definition^^ of an effective local moment, 
(5'?)imp(r) is given by the difference of (5?) of the effective 
SIAM and the medium electron gas without impurity. Since 
the medium changes with temperature, (5'?)imp(7") is calcu- 
lated for each converged medium for the given temperature 
and T* is obtained by the condition (5?)imp(7^*) = 0.05. 

The results for all four different scales are depicted in Fig. 
20 using a fixed filling ntot = 1-6. While the fit to the 
behavior contains a significant error bar due to the logarith- 
mic temperature scale of the NRG and the limited temperature 
range of the Fermi-liquid regime, all other energy scales can 
be obtained with high accuracy due to their definition. The 
dimensionless coupling constant 

2TqU 

has been estimated by the Schrieffer- Wolff transformation.^' 
All four scales show the same exponential dependency on U 
with 7b oc exp(— 1/ag) with a « 0.8 for large U where the 
charge and low-temperature scales are well separated. This 
indicates that, up to a filling-dependent proportionality con- 
stant, there exists only one low-temperature scale. Our find- 
ings are in perfect agreement with the previously reported^^ 
r = estimates of the low-temperature scale but in contrast 
to the Gutzwiller results^^ predicting an enhancement of the 
Tq/Tk ~ exp(l/2g). However, the latter ratio is strongly fill- 
ing dependent as pointed out previously. 

As a consequence of this scaling analysis, {S^){^^{T) is 
collapsed onto one master curve for a fixed filling ntot and 



different values of U using the scaling variable x = T /Tq as 
displayed in Fig. 21. This master curve for (5?)inip(r) shows 
slight deviation from its corresponding shape for the SIAM 
with constant density of states however, it holds up to very 
high temperatures. Clearly visible is the U dependence of the 
maximum of (5?)imp(r). Only in the local moment fixed point 
the value of (5'?)imp = 0.25 is reached. At high temperatures, 
deviations from scaling are related to the crossover from the 
free orbital to the local-moment fixed point^^ and are indepen- 
dent of lattice renormalization effects. 

In addition, we can define a "temperature" dependent 
Kondo scale Tk{T) of the DMFT effective site, by taking the 
DMFT medium density of states Peff(ci)) for a fixed temper- 
ature and solving the condition (5?)imp(7A') = 0.05 for such 
a fixed medium. Since we would find a constant Tk{T) for a 
SIAM, deviations from such a constant function indicate lat- 
tice renormalization effects. This also could give strong hints 
towards the existence of a second low-energy scale. Tk{T 
oo) corresponds to a "free-impurity" Kondo temperature. In 
the Fermi-liquid regime, Tk{T) approaches a constant Tk{Q) 
measuring the lattice low-temperature scale: Tk{0) ~ Tq. In 
Fig. 22, the resulting renormalized dimensionless Kondo tem- 
perature Tk{T)/Tq is plotted versus T /T^ for different values 
of U . Even though two distinguishable plateaus are visible, 
indicating two different low-temperature scales for large val- 
ues of U, they are connected by a fixed ratio independent of 
U. 

The existence of such universal scaling curves for (5'?) imp 
and Tk{T) connecting these different temperature regimes 
provides strong evidence that (i) the low-temperature phase of 
the paramagnetic RAM is characterized by only one thermody- 
namic low-temperature scale, and (ii) all additionally defined 
low-energy scales assigned to features of some physical prop- 
erties are proportional to this single low-temperature scale. 
Our findings are in agreement with the previous DMFT-NRG 
study^^ as well as the series of DMFT-LMA investigations of 
scaling and transport properties. We did not find any in- 
dication for a novel fixed point of the effective site introduced 
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by the nonconstant density of states. In the absence of a phase 
transition, the effective site is characterized by the same fixed 
points as the SIAM^^ and, therefore, only one universal low- 
temperature scale is expected. 

VII. CONCLUSION AND OUTLOOK 

We presented a comprehensive study of the transport prop- 
erties of the periodic Anderson model. We calculated the 
resistivity and the thermoelectric power for metallic heavy- 
fermion systems as well as Kondo insulators using finite tem- 
perature DMFT-NRG and DMFT-MPT approaches. While 
the DMFT-MPT gives extremely reliable results in the weak 
correlation limit p{0)U ^ 1, it fails to reproduce the cor- 
rect low-energy scale for large values of U . In contrast, the 
DMFT-NRG resolves accurately low-energy scales for arbi- 
trary parameters but provides limited spectral resolution in 
finite-temperature spectral functions. The latter has made it 
very challenging to obtain reliable transport properties since 
only excitations of the order T enter the transport integrals. 
Nevertheless, we were able to use DMFT-NRG for the inves- 
tigation of transport properties as well as the four different 
low-energy scales. All these energy scales are related by a 
fixed but filling-dependent ratio, indicating the existence of 
only one relevant low-temperature scale. 

We conclude that all characteristic low-temperature trans- 
port properties of heavy-fermion systems can be explained 
within our numerical DMFT-NRG treatment of the PAM. In 
the metallic situation one obtains a low residual resistivity 
p(0), a p(r) rapidly increasing with increasing temperature 
T following a behavior for very low T <To, a maximum 
of the order 100 fiD. cm at a temperature T^nax (about 10- 
200 K) and a p{T) logarithmically decreasing with increas- 
ing T for T > Tynax, thus showing the typical experimental 
behavior '^''^''^ At the same time the thermoelectric power 
can exhibit several extrema, for some parameters a change 
of sign, and large absolute values of the magnitude of 50- 
150 /iV/K, as seen in experiments. ^^'^^ The optical conduc- 
tivity shows a Drude peak and an additional midinfrared peak 
for low T, whereas these structures merge to a broad structure 
at high T. For Kondo insulators a crossover from the Kondo 
behavior atT > r^ax to an activation behavior is obtained for 
p(r) with a saturation due to the (artificial) finite imaginary 



part 8 (simulating impurity scattering), and the optical con- 
ductivity has no Drude peak for very low T. The absolute 
value of the thermoelectric power S{T) can be even larger in 
or close to the Kondo insulator situation reaching giant values 
of more than 200 /iV/K. The figure of merit can be larger 
than 1, showing that these systems are candidates for thermo- 
electric applications. 

Thus far all our calculations did not include orbital degen- 
eracy which is subject to crystal electric field splittings. These 
additional degrees of freedom are essential for the understand- 
ing of the transport properties of heavy-fermion systems in the 
intermediate temperature regime^^ as well as the broad midin- 
frared peak observed in the optical conductivity.^^ As recently 
pointed out,^^'^^'^'* these additional / configurations can in- 
deed introduce additional low-energy scales due to competi- 
tion or crossover between additional fixed points. Apparently, 
the thermoelectric power is very sensitive to such additional 
configurations since it is determined by particle-hole asym- 
metry around the chemical potential. 

Therefore, in the future this work can and should be ex- 
tended into different directions. On the one hand, the in- 
clusion of realistic /-shell degeneracy and of crystal fields is 
highly desirable to come to a more realistic description of the 
transport properties of heavy-fermion systems valid also for 
high T, where the higher crystal field levels can be excitated. 
On the other hand, the influence of disorder and alloying, in 
particular of the substitution of the (magnetic) lanthanide or 
actinide ions by nonmagnetic ions should be investigated and 
understood theoretically, as there exist many experimental in- 
vestigations systematically studying these substitution effects. 

Note added in proof. After completion of this work, po- 
tentially more accurate ways of calculating the NRG spectral 
functions were proposed^^-^^ which might increase the accu- 
racy of future DMFT-NRG calculations. 
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